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SECTION  1 
INTRODUCTION 


Recent  computational  experience  has  indicated  that  flow 
interactions  with  a  surface  can  materially  alter  precursor  structure 
formation  in  a  non-ideal  airblast.  The  near-surface  precursor  flow 
includes  a  "jet"  that  interacts  with  and  exchanges  momentum  with  the 
surface  boundary  layer.  The  transition  between  the  off-surface  jet 
flow  and  the  flow  in  the  boundary  layer  is  smooth;  the  demarcation 
between  the  two  regions  is  indistinct.  Momentum  lost  to  the  ground 
through  surface  shear  is  diffused  by  turbulence  through  the  boundary 
layer  into  the  jet  without  interruption.  The  jet  velocity  is  reduced 
and  the  jet  pressure  increased  due  to  this  effect.  The  rate  at  which 
work  is  transmitted  to  the  precursor  toe  is  diminished,  resulting  in  a 
reduction  in  the  overall  precursor  length. 

The  detailed  character  of  the  surface  affects  both  the  extent 
of  airflow  momentum  loss  and  the  concomitant  turbulent  diffusivity 
within  the  boundary  flow.  The  collective  surface  conditions  that 
control  these  effects  establish  its  "roughness".  Four  properties  of 
the  surface  appear  to  have  first-order  effects.  These  are:  (1)  the 
type  of  surface  irregularities  (depressions  or  protuberances),  (2)  the 
size  of  the  irregularities,  (3)  the  fractional  amount  of  surface 
covered  by  irregularities,  and  (4)  the  geometry  of  the  irregularities. 

The  effects  of  depression  type  (d-type)  roughness  differ 
somewhat  from  those  of  protuberance-type  (k-type)  roughness  (Reference 
1) .  However,  d-type  roughness  is  considered  to  be  of  lesser  immediate 
interest  and  will  not  be  considered  further  in  this  discussion. 

The  contributions  to  surface  roughness  arise  predominantly  from 
fixed  surface  protrusions.  In  desert  terrains,  fixed  surface 
protrusions  may  include  shrubs,  trees,  rocks,  and  ground  clods. 
Loose,  lightweight  material  such  as  non-fixed  pebbles,  vegetative 
debris,  and  dust  cannot  support  significant  shear.  These  materials 
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tend  to  be  lofted  into  the  boundary  layer  and  diffused  upward  by 
turbulence.  Momentum  is  exchanged  between  lofted  material  and  the 
air,  but  there  is  essentially  no  I oss  of  momentum  to  the  ground,  as 
with  fixed  protrusions. 

A  real  surface  may  comprise  an  essentially  continuous 
distribution  of  protrusion  heights.  Nevertheless,  it  is  useful  to 
discuss  the  isolated  effects  of  groups  cf  these  protrusions  possessing 
grossly  equivalent  geometry  and  size.  In  such  groups,  individual 
members  are  referred  to  as  "roughness  elements".  Typical,  simplified 
roughness  element  models  include  surface-attached  spheres, 
hemispheres,  cones,  and  cylinders. 

Elements  that  are  shorter  than  a  specific  minimum  vertical 
size,  dependent  only  on  flow  conditions,  do  not  contribute  measurably 
to  roughness  effects.  If  a  I  I  roughness  elements  on  a  given  surface 
fall  below  this  height,  the  surface  is  said  to  be  "hydraulically 
smooth" . 

The  steady-state,  clean-flow  (no  dust),  turbulent  boundary 
layer  of  a  flow  over  a  hydraulically  smooth  surface  consists  of  four 
distinctive  regions.  The  innermost  of  these  regions  contacts  the 
surface  and  extends  vertically  to  about  0.1%  of  the  overall  boundary- 
layer  thickness.  The  flow  in  this  "sublayer"  is  essentially  laminar 
and  the  velocity  increases  linearly  with  height  from  zero  at  the 
surface  Turbulence  kinetic  energy  grows  rapidly  above  the  laminar 
sublayer  in  a  second  region  referred  to  as  the  "buffer"  layer.  This 
growth  levels  off  at  about  1-2%  of  the  boundary  layer  thickness,  which 
marks  the  beginning  of  the  third,  "log-law"  layer.  The  turbulence 
energy  is  nearly  constant  throughout  the  log- 1  aw  region  which  extends 
to  about  10%  of  the  boundary  layer,  thickness.  Flow  conditions  in  the 
log-law  region  are  well  behaved.  Turbulence  production  and 
dissipation  are  in  a  state  of  near  equilibrium  which  is  reflected  by  a 
relatively  slow  variation  in  turbulence  kinetic  energy.  The  velocity 
increases  logarithmically  with  height. 
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The  iaminar  sublayer  and  buffer  layer  taken  together  are  often 
referred  to  as  the  "viscous  sublayer".  The  viscous  sublayer  and  the 
log- 1  aw  region  together  comprise  the  "inner  layer"  of  the  boundary 
layer  flow.  The  remaining  90fl  of  the  boundary  layer  comprises  the 
"outer  region".  In  simple  flat-plate  flow,  turbulence  energy  drops 
off  in  the  outer  region  while  the  turbulence  length  scale  increases. 

The  four-layer  smooth-wall  boundary  flow  structure  begins  to  be 
disrupted  when  surface  roughness  elements  protrude  above  the  laminar 
sublayer  into  the  buffer  region.  The  extent  of  this  disruption 
increases  with  further  increases  in  roughness  element  size. 
Turbulence  near  the  surface  is  enhanced  and  small  turbulent  eddies 
penetrate  into  the  laminar  sublayer.  When  the  roughness  elements  are 
large  enough  as  to  extend  into  the  log- law  region,  a  viscous  sublayer 
is  completely  absent.  The  inner  layer  then  consists  only  of  a  log- 1  aw 
region,  and  the  surface  is  said  to  be  "fully  rough". 

A  criterion  widely  used  to  distinguish  a  hydraulically  smooth 
surface  is 


where  yp  is  the  roughness  element  height,  v  is  the  kinematic  viscosity 
of  the  flow,  and  uT  is  the  "shear"  or  "friction"  velocity,  uT  = 
(rw/Pw)1/2  (,n  *hich  rw  is  the  surface  shear  stress  and  is  the 
fluid  density  at  the  wall).  On  the  other  hand,  the  surface  can  be 
considered  fully  rough  when 

yR>70^.  (2) 

These  criteria  cannot  easily  be  applied,  unfortunately,  unless 
an  estimate  (or  measurement)  of  the  surface  shear  is  available.  For 
some  purposes,  a  satisfactory  estimate  for  this  purpose  can  be  taken 
by  reference  to  an  approximately  comparable  flat  plate  flow.  In  a 
simple  incompressible  flow  over  a  flat  plate,  the  skin  friction 
coefficient,  Cf  =  Tw/(0.5pu^)  can  be  found  from  the  Schoenherr 
formula, 
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Cf  1/2[l -80  £ji(RexCf)  ♦  1.7]  =  1  , 


(3) 


where  Rex  is  the  Reynolds  number  based  on  the  distance  from  the 
leading  edge  of  the  plate,  Rex  =  u^X/t/,  and  is  the  velocity 

external  to  the  boundary  layer.  A  useful  direct  relation  between  ur 
and  Cf  is 

—  =  (C,/2)1/2  (4) 

U^  ^  f 
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SECTION  2 

EFFECT  OF  AREAL  COVERAGE 

The  remainder  of  the  present  discussion  will  consider  only  the 
fully  rough  condition.  For  this  case,  the  velocity  profile  of  a 
simple,  incompressible,  non-separating  boundary  layer,  up  to  about  1055 
of  the  boundary  layer  thickness,  is  well  represented  by  the  rough-wall 
Clauser  "wall  function". 


u 

u 


T 


♦  C' 


(5) 


Here  x  =  0.41  (Von  Karman’s  constant)  is  fixed  and  universal,  but  the 
value  of  C'  varies  with  the  roughness  element  geometry  and  area 
coverage . 


Figure  1  is  a  plot  of  C' 
coverage, 


as  a  function  of  fractional  area 


fR  =  Ar/A  , 


(6) 


where  Ar  is  the  horizontal  projected  surface  area  covered  by  roughness 
within  area  A.  Data  points  on  this  figure  were  derived  from  Reference 
2.  The  data  points  shown  individually  on  this  plot  are  for  three- 
dimensional  roughness  elements  as  opposed  to  two-dimensional  elements 
such  as  transverse  rods. 

Smaller  values  of  C'  reflect  greater  ground  shear  on  surfaces 
having  roughness  elements  of  equal  height  yR.  This  is  true  regardless 
of  changes  in  the  geometry  of  the  elements,  if  the  exterior  flowfield 
is  held  fixed.  The  data  for  spherical  roughness  elements  show  that  C' 
is  not  a  monotonic  function  of  f r .  A  minimum  occurs  at  20-25% 
(depending  on  which  curvefit  is  used).  The  data  for  other  three- 
dimensional  roughness  elements  is  insufficient  on  Figure  1  to  exhibit 
a  minimum.  However,  the  curvefit  by  Dvorak  (Reference  3),  which  was 
fit  through  data  for  both  two-  and  three-dimensional  roughness 


5 


auser  Log-Law  constant  C  vs.  fractional  roughness  coverage,  fp  =  AR/A. 
R  =  area  covered  by  roughness,  A  =  total  area;  XR  =  l/fR  is  "roughness 


elements,  and  that  by  Betterman  (Reference  4)  for  strictly  two- 
dimensional  elements  are  also  shown  on  the  figure.  The  collective 
data  to  which  these  curves  were  fit  suggest  that  a  minimum  value  of  C' 
exists  for  all  roughness  geometries,  but  that  the  value  of  the  minimum 
may  change  according  to  the  geometry. 

The  drag  on  a  single  roughness  element  is  a  function  of  its 
vertical  projected  area,  or  "lateral  area",  A^,  while  the  amount  of 
surface  covered  by  one  element  is  a  function  of  its  hor i zonta I 
projected  area,  Ap ' .  For  example,  a  vertical  cylinder  of  diameter  d 
and  height  h  has  a  lateral  area  of  A|_  =  dh  and  a  specific  surface  area 
of  Ap'  =  *dz/4 .  It  is  not  surpr  ising  that  the  values  of  C'  are 
smaller  for  those  roughness  element  geometries  which  have  larger 
values  of  A(_/Ap'.  For  the  surface-tangent  spheres  on  Figure  1,  A^/Ap' 
=  1.  By  comparison,  the  conical  elements  (yp/d  =  0.469)  have  A|_/Ar'  = 
0.299,  and  the  angles  (neglecting  their  horizontal  part,  and  assuming 
the  thickness  of  the  vertical  part  to  be  d/5),  have  A^/Ar'  =  1.88. 

As  the  number  of  roughness  elements  per  unit  area  (that  is,  the 
fractional  coverage)  increases,  the  drag  coefficient  associated  with 
the  individual  elements  initially  remains  constant.  However,  when  the 
spacing  between  elements  is  reduced  to  several  element  heights,  flow 
interactions  between  elements  cause  the  velocity  in  their  vicinity  to 
decrease.  The  drag  coefficient  based  on  the  exterior  flow  velocity 
then  abruptly  falls  off.  This  effect  has  been  investigated  in  some 
detail  by  Marshall  (Reference  5),  and  helps  to  explain  the 
characteristic  shape  of  the  C'  vs.  fR  curves  in  Figure  1.  When  fR  is 
small,  the  drag  coefficient  of  individual  elements  is  at  its  highest 
level  and  does  not  change  appreciably  with  moderate  increases  in  fp. 
Hence  a  larger  number  of  roughness  elements  simply  transmits  a  larger 
total  force  to  the  surface.  This  is  reflected  by  a  steady  decrease  in 
C'  with  increasing  fp.  Above  the  coverage  at  which  the  element  drag 
coefficient  begins  to  decline,  increases  in  the  number  of  elements  is 
offset  by  decreases  in  the  drag  per  element.  The  total  force 
transmitted  to  the  surface  drops  off  with  further  increases  in  fp,  and 
the  value  of  C'  increases. 


It  is  of  interest  to  note  that  in  Marshall’s  experiments,  only 
the  amount  of  coverage  and  not  the  d i str i but i on  of  roughness  elements 
on  the  surface  appeared  to  influence  the  results.  That  is,  regular 
arrays  and  randomly  distributed  elements  produced  equivalent  effects 
in  these  experiments.  Most  of  the  tests  were  carried  out  with 
cylindrical  roughness  elements  of  varying  Aj_/Ap '  (varying  d/h) , 
although  some  experiments  were  conducted  with  hemispheres.  Another 
qualification  to  note  is  that  Marshall’s  tests,  in  common  with  all 
other  experiments  referenced  in  this  paper, were  carried  out  in  the  low 
subsonic  regime. 
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SECTION  3 

EQUIVALENT  SURFACE  ROUGHNESS 


Due  to  the  character  of  the  C'  vs.  fp  (and  geometry)  curves,  the 

relative  roughness  effects  of  different  surfaces  cannot  be  compared  in 

terms  of  a  single,  directly  measurable  quantity.  A  classical  basis 
for  comparison  is  in  terms  of  "equivalent  sand-grain  roughness  size", 
Ks .  This  is  the  diameter  of  fixed  sand- 1  ike  elements,  packed  on  a 

surface  such  that  fps  =  0.91  (and  C'  =  8.5),  which  would  yield  the 

same  wall  shear  as  that  on  a  specimen  surface  with  a  particular 
roughness  element  geometry,  height,  and  coverage.  If  the  value  of  Ks 
is  known  for  the  specimen  surface,  the  Clauser  wall  function,  Eqn.(5), 
would  be  written  as 

r  *  i  «"GH  * 8  5  •  <?> 

T  S 

In  this  way,  Ks  becomes  a  measure  of  "roughness  scale". 

An  alternative,  equivalent  procedure  is  to  rewrite  Eqn.  (5)  in 
the  form 


where  roughness  scale  Z0  is  related  to  C'  and  yp  through 

ZQ  =  yp  EXP  (-*C')  ■  (9) 

Both  Eqn.  (7)  and  Eqn.  (8)  appear  in  the  literature. 

Surfaces  with  equal  values  of  Ks  or  with  equal  values  of  Z0  are 
hydrodynamica I ly  equivalent  in  terms  of  roughness  effects.  Choosing 
Eqn.  (8)  as  the  basis  for  present  comparisons,  the  Dvorak  fit,  along 
with  the  two-point  fit  on  Figure  1  (for  surface-tangent  spheres,  above 
fp  =  0.247),  becomes 

(1.36  ypfR117  ,  fR  <  0.247 

zo  =  R  R  R  (10) 

.0.0158  yRfR'201  ,  fp  >  0.247 
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For  example,  the  value  of  C'  for  a  surface  covered  by  2  cm.  dia.  rocks 
j  (yp  =  2  cm)  over  218  of  the  surface  area  (fp  =  0.21),  is  C'  =3.65. 

The  value  of  C'  for  a  terrain  covered  by  100  cm  dia.,  dense  spherical 
shrubs  (yp  =  100  cm)  over  0.748  of  the  surface  area  (fp  =  0.0074)  is 
C'  =  6.67.  Both  of  these  conditions  correspond  to  ZQ  =  0.438  cm  in 
|  Eqn.  (8).  In  terms  of  hydrodynamic  roughness  effects,  the  two  surfaces 

are  equivalent. 

The  example  cited  is  significant  because  it  illustrates  that  a 
terrain  in  which  the  surface  appears  to  be  relatively  smooth  but  is 

vegetated  with  small,  dense  shrubs  scattered  over  about  18  of  the 

ground  area  is  hydrodynamics  I ly  as  rough  as  if  the  surface  were 

"heavily  covered"  with  rocks.  The  separation  between  2  cm  dia 
spherical  "rocks",  corresponding  to  an  areal  surface  coverage  of  218, 
is  only  about  4.2  cm  in  a  uniform  hexagonal  array. 

In  order  to  characterize  a  surface  for  boundary  layer 

calculations,  only  the  value  of  ZQ  (or  Ks)  is  needed.  However,  Eqn 
(10)  applies  approximately  only  for  spherical  roughness  elements  and 
then  only  at  relatively  low  Mach  and  Reynolds  numbers.  It  appears 
that  the  effects  of  higher  flow  speeds  on  the  values  of  ZQ  have  not 
been  investigated  experimentally.  Flow  regime  effects  seem  almost 
certain  to  influence  ZQ  through  the  impact  of  Mach  and  Reynolds 
numbers  on  the  drag  coefficients  of  individual  roughness  elements. 
Without  detailed  experimental  or  calculated  information  of  this  kind, 
it  is  difficult  to  estimate  appropriate  values  for  Z0  from  field 
measurements,  in  a  given  terrain,  of  yp,  fp  and  roughness  element 
geometry . 
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SECTION  4 

DISTRIBUTED  ROUGHNESS  SIZES 


An  additional  complication  arises  when  roughness  elements  of 
different  sizes  and  geometries  are  present.  This  is,  of  course,  the 
case  for  all  real  surfaces,  although  the  effects  of  one  group  of 
roughness  elements  may  dominate.  As  an  example,  a  hypothetical  desert 
surface  is  offered  in  Table  1.  On  this  hypothetical  surface,  loose 
dust  has  not  been  included.  This  is  because  the  magnitude  of  shear 
that  loose  dust  can  sustain  (before  its  movement  is  initiated)  is  very 
small  compared  to  the  shear  transmitted  by  larger,  fixed  roughness 
elements.  The  hypothetical  surface  includes  roughness  elements 
spanning  four  orders  of  magnitude.  All  are  assumed  to  be  spherical. 
The  majority  of  the  surface  (805t  by  area)  is  covered  by  fixed,  small 
"sand  grains*.  Nevertheless,  the  major  contributors  to  roughness  are 
"medium  rocks"  and  "small  shrubs",  each  of  which  covers  only  IS  of  the 
surface.  It  would  probably  be  satisfactory  for  engineering 
calculations  to  use  ZQ  =  0.622  cm.  in  Eqn .  (8)  for  this  surface  (or 
perhaps  a  slightly  smaller  value). 

It  is  worth  emphasizing  that  ZG  bears  little  direct  relation  to 
the  physical  size  of  roughness  elements.  The  actual  height  of  the 
elements  is  yp,  while  ZQ  is  a  characteristic  roughness  scale  that 
accounts  for  element  size,  geometry,  and  areal  coverage. 

It  may  be  possible  to  establish  the  value  of  ZQ  for  a  given 
terrain  without  a  detailed  survey  like  that  given  in  Table  1.  An 
example  of  this  appears  in  Reference  6.  The  velocity  profile  of  an 
atmospheric  wind- induced  boundary  layer  was  measured  (at  the  GMX  area 
of  the  Nevada  Test  Site).  The  measurements  were  made  using  a  pitot- 
static  rake  of  appropriate  height.  From  the  velocity  profile  it  was 
possible  to  establish  a  characteristic  roughness  scale. 


Table  1. 


Distribution  of  roughness  scale  Z0  for  hypothetical 
desert  surface.  Values  for  ZQ  assume  spherical  roughness 
elements  of  diameter  yg  and  fractional  areal  coverage  fp. 


ROUGHNESS  ELEMENT  TYPE 

>R 

(cm) 

fR 

— r0 - 

(cm) 

Sma 1 1  sand  gra i ns 

0.02 

0.80 

0.000495 

Large  sand  grains 

0.2 

0.16 

0.0319 

Sma 1 1  rocks 

2.0 

0.02 

0.0280 

Medium  rocks 

10. 

0.01 

0.0622 

Small,  dense  shrubs 

100. 

0.01 

0.622 

—i 

1 


i 

4 


_ i 


While  this  straightforward  procedure  is  relatively  convenient, 
it  is  certainly  not  clear  that  the  roughness  scale  determined  from 
such  measurements  can  be  applied  unmodified  to  the  boundary  layer  of  a 
high-speed,  transient  flow.  It  would  accordingly  seem  worthwhile  to 
investigate  the  possibility  for  establishing  procedures  from  which  the 
correct  scale  associated  with  a  high-speed  flow  could  be  inferred  from 
the  value  determined  through  low-speed  wind  measurements. 


As  an  alternative,  it  may  prove  feasible  to  establish  ZQ  for  a 
specific  terrain  directly  from  boundary  layer  profiles  measured  during 
an  airblast  event.  The  smallest  scale  at  which  measurements  of  this 
kind  would  be  meaningful  is  uncertain,  but  an  estimate  might  be  RQ  » 
yp,  where  RQ  is  a  characteristic  airblast  radius  for  a  surface  burst, 


R  = 
o 


(11) 


in  which  EQ  is  the  explosion  energy  and  Pj  is  initial  ambient 
pressure.*  For  example,  with  RQ  =  20yp,  Pj  =  0.1  MPa,  and  yp  =  1  m, 
the  smallest  explosion  energy  would  be  EQ  =  5.02  x  10®  J,  This 
suggests  that  a  0.0024  KT  (or  2.4T)  surface  burst  might  reasonably 
suffice  for  a  measurement  of  ZQ. 


*  At  the  ground  range  for  which  R  =  RQ,  the  overpressure  is  AP  =  1.4 
P1 
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SECTION  5 
ROUGHNESS  SCALING 


Finally,  it  is  noted  that  the  characteristic  roughness  Z0  is  a 
fixed  property  of  the  terrain  and,  perhaps,  of  the  flow  regime  (Mach 
and  Reynolds  numbers).  As  such,  the  value  of  ZQ  for  a  particular 
surface  does  not  change  as  the  airbfast  scale  is  changed;  i.e.,  in 
this  sense  ZQ  does  not  "scale".  However,  for  some  assessments  of 
airblasts  at  different  scales,  it  is  of  interest  to  hold  the  global 
effects  of  surface  roughness  constant.  In  order  for  the  overall 
effects  of  roughness  to  remain  unchanged  as  the  airblast  scale  is 
changed,  it  is  necessary  to  allow  Z0  to  vary.  Experimentally,  this 
would  require  tests  to  be  conducted  over  surfaces  with  different  Z0  as 
E0  is  changed.  It  can  be  shown  that  the  appropriate  scaling  for  equal 
roughness  effects  is  Z^/Z^  =  (E02/E0i)*^.  (That  is,  "cube-root" 
scaling  is  valid.)  For  example,  to  simulate  the  roughness  effects  of 
the  hypothetical  desert  surface  in  Table  1  at  1:27,000  yield  scale 
(1:30  length  scale),  155  of  the  test  surface  should  be  covered  with 
fixed  spheres,  3.3  cm.  in  diameter.  Experiments  at  this  scale  are 
within  the  capability  of  existing  confined-flow  airblast  simulators. 


14 


SECTION  6 

LIST  OF  REFERENCES 

1.  Perry, A.  E.,  et.  al.,  "Rough  Wall  Turbulent  Boundary  Layers,"  J. 
Fluid  Mech.  Vol .  37,  Part  2,  1969,  pp .  383-413. 

2.  Schlichting,  H.,  Boundary  Layer  Theory,  McGraw-Hill,  New  York, 
1968. 

3.  Betterman,  D.,  "Contribution  A  L’ Etude  de  la  Conche  Limite 
Turbulante  Le  Long  De  Plaques  Regueuses,"  Rapport  65-6,  Centre 
Nat.  de  la  Rech.  Sc i . ,  Paris  France,  1965.  Referenced  in  Cebeci, 
T.  and  Smith,  A.M.O.,  Analysis  of  Turbulent  Boundary  Layers,  New 
York.  Academic  Press,  1974,  p.  130. 

4.  Dvorak,  F.A.,  "Calculation  of  Turbulent  Boundary  Layers  on  Rough 
Surfaces  in  Pressure  Gradient,"  AIAA  J.,  Vol.  7,  No.  9,  1969, 

p.  1752. 

5.  Marshall,  J.K.,  "Drag  Measurements  in  Roughness  Arrays  of 
Varying  Density  and  Distribution",  Agr.  Meteorol.,  Vol.  8,  1971, 
pp.  269-292. 

6.  Shinn,  J.  H.,  et  a  I ,  "Observations  of  Dust  Flux  in  the  Surface 
Boundary  Layer  for  Steady  and  Nonsteady  Cases,"  in  Atmosphere 
Surface  Exchange  of  Particulate  and  Gaseous  Pollutants  ^1974^7 
Englemann,  R.J  ,  and  Sehmel,  G.A.,  Eds.,  NTIS  CONFERENCE-740921, 
January.  1976,  p.  625. 


15 


DISTRIBUTION  LIST 

DNA-TR-88-140 


DEPARTMENT  OF  DEFENSE 

AFSOUTH 

ATTN:  U  S  DOCUMENTS  OFFICER 

ARMED  FORCES  STAFF  COLLEGE 
ATTN.  LIBRARY 

ASSISTANT  TO  THE  SECRETARY  OF  DEFENSE 

ATOMIC  ENERGY 

ATTN:  EXECUTIVE  ASSISTANT 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  DB 
ATTN:  G FERRELL 
ATTN:  R  MANN 
ATTN:  DE(ESTIMATES) 

ATTN:  DI-5 

ATTN:  DIA/VPA-2  (FED  RES  DIV) 

2  CYS  ATTN:  RTS-2B 

ATTN:  RTS-2C  (TECH  SVCS  &  SPT) 

DEFENSE  NUCLEAR  AGENCY 
ATTN:  DFRA 
ATTN:  NANF 
ATTN:  NASF 
ATTN:  NAWE 
ATTN:  OPNA 
ATTN:  OPNS 
ATTN  RAEE 
4  CYS  ATTN:  TITL 

DEFENSE  NUCLEAR  AGENCY 
ATTN:  NVCG 

DEFENSE  NUCLEAR  AGENCY 

ATTN:  TDTTWSUMMA 

DEFENSE  TECHNICAL  INFO  CENTER 
2  CYS  ATTN:  DTIC/FDAB 

DIRECTOR  NET  ASSESSMENTS 

ATTN:  DOCUMENT  CONTROL 

INTELLIGENCE  CENTER,  PACIFIC 
ATTN:  COMIPAC 

JOINT  DATA  SYSTEM  SUPPORT  CTR 
ATTN:  R  MASON 
ATTN:  C-332 

JOINT  STRAT  TGT  PLANNING  STAFF 
2  CYS  ATTN:  JK  (ATTN:  DNA  REP) 

ATTN:  STUKMILLER 
ATTN:  JLT 
ATTN:  JP 
ATTN:  JPEP 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  DNA-LL 

NATIONAL  DEFENSE  UNIVERSITY 
ATTN.  ICAF  ICC 
ATTN:  NWCO 


NATIONAL  SECURITY  AGENCY 
ATTN:  CHIEF  A  GROUP 

OFFICE  OF  THE  JOINT  CHIEFS  OF  STAFF 
ATTN:  ED30(J-3  STRAT  OPS  DIV) 

ATTN:  J-3/NUC  OPNS  BR,  STRAT  OPNS  BR 
ATTN:  J-5  NUCLEAR  &  CHEM  DIV 
3  CYS  ATTN:  J-5  NUC  DIV/STRAT  DIV/FP&P  DIV 
ATTN:  J-8/CAD 
ATTN.  JAD/SFD 
ATTN:  JAD/SSD 

OFFICE  OF  THE  SEC  OF  DEFENSE 
ATTN:  NAVAL  FORCES 
2  CYS  ATTN:  STRAT  PROGRAMS  &  TNF 

STRATEGIC  AND  THEATER  NUCLEAR  FORCES 
ATTN:  DR  WOODRUFF 

USEURC  EAN  COMMAND 
ATTN:  ECC3S-CCN 

U  S  EUROPEAN  COMMAND 
ATTN:  ECJ-LW 

U  S  EUROPEAN  COMMAND 
ATTN:  ECJ-2-ITD 

U  S  EUROPEAN  COMMAND 
ATTN:  ECJ-3 

U  S  EUROPEAN  COMMAND 
ATTN:  ECJ2-T 

U  S  EUROPEAN  COMMAND 
ATTN:  ECJ5-N 

U  S  NATIONAL  MILITARY  REPRESENTATIVE 
ATTN:  U  S  DOCUMENTS  OFFICER 

UNDER  SEC  OF  DEFENSE  FOR  POLICY 
ATTN:  DUSP/P 
ATTN:  USD/P 

UNDER  SECRETARY  OF  DEFENSE 

ATTN.  DEP  UND  SEC,  TAC  WARFARE  PROG 
ATTN:  GSEVIN 

DEPARTMENT  OF  THE  ARMY 

DEP  CH  OF  STAFF  FOR  OPS  &  PLANS 
ATTN.  DAMO-NCN 
ATTN:  DAMO-NCN 
ATTN:  DAMO-RQS 
ATTN:  DAMO-ZXA 

DEP  CH  OF  STAFF  FOR  RSCH  DEV  &  ACQ 
ATTN:  DAMA-CSM-N 

EIGHTH  U  S  ARMY 

ATTN:  EACJ-PON-NS 

HARRY  DIAMOND  LABORATORIES 
ATTN:  SLCIS-IM-TL 


DNA-TR-88-140  (DL  CONTINUED) 


U  S  ARMY  AIR  DEFENSE  SCHOOL 
ATTN:  COMMANDANT 

U  S  ARMY  ARMOR  SCHOOL 
ATTN:  ATSB  CTD 
ATTN:  TECH  LIBRARY 

U  S  ARMY  BALLISTIC  RESEARCH  LAB 
ATTN:  SLCBR  DD  T 
ATTN:  SLCBR  SS  T 

U  S  ARMY  COMB  ARMS  COMBAT  DEV  ACTY 
ATTN:  ATZL  CAP 

U  S  ARMY  COMD  &  GENERAL  STAFF  COL 
ATTN:  ACQ  LIBRARY  DIV 
ATTN:  ATZL  SWT  A 

U  S  ARMY  EUROPE  AND  SEVENTH  ARMY 
ATTN:  AEAGC 

U  S  ARMY  INFANTRY  CTR  &  SCH 
ATTN:  ATSH-CD-CS 

U  S  ARMY  MATERIEL  COMMAND 
ATTN  AMCCN 

U  S  ARMY  MATERIEL  SYS  ANALYSIS  ACTVY 
ATTN  AMXSYCR 

U  S  ARMY  MISSILE  COMMAND 
ATTN:  AMSMI  XF 

U  S  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
ATTN  MONA  NU 

U  S  ARMY  NUCLEAR  EFFECTS  LABORATORY 
ATTN:  DRJMEASON 

U  S  ARMY  WAR  COLLEGE 
ATTN:  LIBRARY 

USA  SURVIVABILITY  MANAGMENT  OFFICE 
ATTN:  J BRAND 

DEPARTMENT  OF  THE  NAVY 

MARINE  CORPS  DEV  &  EDUCATION  COMMAND 
ATTN:  COMMANDER 

NAVAL  AIR  FORCE 

ATTN:  COMMANDER 

NAVAL  AIR  FORCE 


NAVAL  RESEARCH  LABORATORY 
ATTN:  CODE  1220 
ATTN:  CODE  2627  (TECH  LIB) 

NAVAL  SEA  SYSTEMS  COMMAND 
ATTN:  SEA-09G53  (LIB) 

ATTN:  SEA-643 

NAVAL  SURFACE  FORCE 

ATTN:  COMMANDER 

NAVAL  SURFACE  FORCE 

ATTN:  COMMANDER 

NAVAL  WEAPONS  EVALUATION  FACILITY 
ATTN:  CLASSIFIED  LIBRARY 

NUCLEAR  WEAPONS  TNG  GROUP,  ATLANTIC 
ATTN:  CODE  221 
ATTN:  DOCUMENT  CONTROL 

NUCLEAR  WEAPONS  TNG  GROUP,  PACIFIC 
ATTN,  DOCUMENT  CONTROL 

OFC  OF  THE  DEP  CHIEF  OF  NAVAL  OPS 
ATTN.  OP  654 
ATTN:  OP  981 

OFFICE  OF  THE  CHIEF  OF  NAVAL  OPNS 
ATTN:  CNO  EXECUTIVE  PANEL 

OPERATIONAL  TEST  &  EVALUATION  FORCE 
ATTN:  CODE  80 

OPERATIONAL  TEST  &  EVALUATION  FORCE, 
ATTN:  INTEL  OFFICER 

TACTICAL  TRAINING  GROUP,  PACIFIC 
ATTN:  COMMANDER 

TACTICAL  WINGS  ATLANTIC 
ATTN:  COMMANDER 

THEATER  NUCLEAR  WARFARE  PROGRAM  OFC 
ATTN:  PMS423 

DEPARTMENT  OF  THE  AIR  FORCE 

AFIA/INKD 

ATTN  MAJ  COOK 

AFIS/INT 

ATTN:  INT 


ATTN:  COMMANDER  AIR  FORCE  SYSTEMS  COMMAND 

ATTN:  DL 

NAVAL  INTELLIGENCE  SUPPORT  CTR  ATTN  SD 

ATTN:  NISC  30 

AIR  FORCE  WEAPONS  LABORATORY 

NAVAL  OCEAN  SYSTEMS  CENTER  ATTN  SUL 

ATTN:  CODE  9642  (TECH  LIB) 

AIR  UNIVERSITY  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL  ATTN  AUL-LSE 

ATTN:  CODE  1424  LIBRARY 

ASSISTANT  CHIEF  OF  THE  AIR  FORCE 
ATTN:  SAF/ALR 


Dist-2 


DNA-TR-88-140  (OL  CONTINUED) 


DEPUTY  CHIEF  OF  STAFF/AF-RDQM 
ATTN.  AF/RDQI 

DEPUTY  CHIEF  OF  STAFF/XOO 
ATTN:  AF/XOC 

DEPUTY  CHIEF  OF  STAFF/XOX 
ATTN:  AFXOXFM 
ATTN.  AFXOXFS 

SPACE  DIVISION/CWH 
ATTN:  CWH 

STRATEGIC  AIR  COMMAND 
ATTN:  NRI/STINFO 

STRATEGIC  AIR  COMMAND 

ATTN:  STIC  (544SIW) 

TACTICAL  AIR  COMMAND 
ATTN:  TAC/DOA 

TACTICAL  AIR  COMMAND 
ATTN:  TAC/XPJ 

U  S  AIR  FORCES  IN  EUROPE/DEX 
ATTN:  USAFE/DEXX 

U  S  AIR  FORCES  IN  EUROPE/DOT 
ATTN:  USAFE/DOQ 

US  AIR  FORCES  IN  EUROPE 
ATTN:  USAFE/INAT 

DEPARTMENT  OF  ENERGY 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  L-35 
ATTN:  L  38 
ATTN:  L-389 
ATTN:  W HOGAN 
ATTN:  L-53  TECH  INFO  DEPT  LIB 
ATTN:  Z  DIVISION  LIBRARY 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  D  RICHMOND 
ATTN:  TDOWLER 
ATTN:  REPORT  LIBRARY 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  TECH  LIB  3141 
ATTN:  RB  STRATTON 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  NIO 

FEDERAL  EMERGENCY  MANAGEMENT  AGENCY 
ATTN.  GORRELL 
ATTN:  JF  JACOBS 

U  S  DEPARTMENT  OF  STATE 
ATTN:  PM/TMP 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ADVANCED  RESEARCH  &  APPLICATIONS  CORP 
ATTN:  DOCUMENT  CONTROL 

RDM  CORP 

ATTN.  CWASAFF 
ATTN:  J  BODE 
ATTN:  J  BRADDOCK 
ATTN.  R BUCHANAN 

BOEING  CO 

ATTN:  H  WICKLEIN 
ATTN:  JW  RUSSELL 

COMPUTER  SCIENCES  CORPORATION 
ATTN:  F  EISENBARTH 

GRUMMAN  DATA  SYSTEMS  CORPORATION 
ATTN:  SSHRIER 

HORIZONS  TECHNOLOGY,  INC 
ATTN:  J  PALMER 

I  IT  RESEARCH  INSTITUTE 

ATTN:  DOCUMENTS  LIBRARY 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  CLASSIFIED  LIBRARY 
ATTN:  JGROTE 

JAYCOR 

ATTN:  R  SULLIVAN 

KAMAN  SCIENCES  CORP 
ATTN:  F  SHELTON 

KAMAN  SCIENCES  CORP 
ATTN:  E  CONRAD 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

KAMAN  TEMPO 

ATTN:  DASIAC 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN.  J  DONATHAN 

ORION  RESEARCH  INC 

ATTN:  JESCHOLZ 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  HBRODE 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  D  GORMLEY 

R  &  D  ASSOCIATES 

ATTN:  C  K  B  LEE 
ATTN:  D  SIMONS 

2  CYS  ATTN:  DOCUMENT  CONTROL 

R  &  D  ASSOCIATES 

ATTN:  C  KNOWLES 
ATTN:  J  THOMPSON 


Dist-3 


DNATR-88-140  (DL  CONTINUED) 


R  &  D  ASSOCIATES 

ATTN:  GGANONG 

RAND  CORP 

ATTN:  V  JACKSON 

RAND  CORP 

ATTN:  B BENNETT 

SCUBED 

ATTN:  BPYATT 
2  CYS  ATTN:  T  PIERCE 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DOCUMENT  CONTROL 
ATTN:  ESWICK 
ATTN:  J  MARTIN 
ATTN:  M  DRAKE 
ATTN:  R  J  BEYSTER 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  SHANNON 
ATTN:  LGOURE 
ATTN:  WLAYSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  RCRAVER 

TRW  INC 

ATTN:  DSCALLY 
ATTN:  R  BURNETT 


Dist-4 


